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Recent experiments on the anisotropic spin- 1/2 triangular antiferromagnet Cs2CuBr4 have re- 
vealed a remarkably rich phase diagram in applied magnetic fields, consisting of an unexpectedly 
large number of ordered phases. Motivated by this finding, we study the role of three ingredients — 
spatial anisotropy, Dzyaloshinskii-Moriya interactions, and quantum fluctuations — on the magneti- 
zation process of a triangular antiferromagnet, coming from the semiclassical limit. The richness of 
the problem stems from two key facts: 1) the classical isotropic model with a magnetic field exhibits 
a large accidental ground state degeneracy, and 2) these three ingredients compete with one another 
and split this degeneracy in opposing ways. Using a variety of complementary approaches, including 
extensive Monte Carlo numerics, spin-wave theory, and an analysis of Bose-Einstein condensation of 
magnons at high fields, we find that their interplay gives rise to a complex phase diagram consisting of 
numerous incommensurate and commensurate phases. Our results shed light on the observed phase 
diagram for Cs2CuBr4 and suggest a number of future theoretical and experimental directions that 
will be useful for obtaining a complete understanding of this material's interesting phenomenology. 
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I. INTRODUCTION 

Antiferromagnetic spin models on the triangular lattice 
constitute one of the simplest and most widely studied 
realizations of geometric frustration. Indeed, the Ising 
triangular antiferromagnet was the first spin model found 
to possess a disordered ground state and extensive resid- 
ual entropyi- at zero temperature (T) . While the classical 
Heisenberg model on the triangular lattice does order at 
T = into a well-known commensurate spiral pattern 
(also known as a a/S x \fZ state), the fate of the quan- 
tum spin-1/2 Heisenberg Hamiltonian has been the sub- 
ject of a long and fruitful debate spanning over 30 years 
of research. Although the originally proposed resonat- 
ing valence bond liquic? ^ did not emerge as the ground 
state of the spin-1/2 Heisenberg modeP"^!, such a phase 
was later found in a related quantum dinier model on the 
triangular lattice^. 

Triangular antiferromagnets in an applied magnetic 
field — which is our focus here — have also been extensively 
studied for decades, and found to possess unusual mag- 
netization physics which remains only partially under- 
stood. Underlying much of this interesting behavior is 
the discovery, made long agcP, that in a magnetic field 
Heisenberg spins with isotropic exchange interactions ex- 
hibit a large accidental classical ground state degeneracy. 
That is, at finite magnetic fields there exists an infinite 
number of continuously deformable classical spin config- 
urations which constitute minimum energy states, but 
are in no way symmetry related. As reviewed later, this 
degeneracy is lifted by thermal (finite Tpl and quantum 
(finite spin S y fluctuations; such fiuct uat io n- driven se- 
lection, also known as order-by-disordep^^I^, results in a 
nontrivial temperature-field phase diagrarrP consisting of 
three ordered phases. In particular, coplanar 'Y' and 'V 



states are separated by a coUinear up-up-down (UUD) 
phase which realizes, for finite spin S and at T = 0, a 
one-third magnetization plateau over a finite field inter- 
val. This plateau state preserves continuous spin rota- 
tion symmetry about the magnetic field direction and is 
remarkably stable: unlike all other magnetically-ordered 
states it survives "dimensional reduction" and exists even 
in the smallest pos sible t riangular lattice strip, the two- 
chain zig-zag laddepi^ESl, 

To date a large number of magnetic materials have 
been synthesized that realize triangular antiferromag- 
nets, and experiments on such compounds have high- 
lighted the spectacular breadth of phenomena that can 
be driven by the interplay between magnetic fields and 
geometric frustration stemming from the lattice. For ex- 
ample, a stacked triangular antiferromagnet with weak 
inter-plane coupling is realized by the S = 5/2 ma- 
terial RbFe(Mo04)2, whose phase diagrarrfSl features 
all of the fluctuation-selected states discussed above. 
Reducing the magnitude of the magnetic ion's spin 
enhances quantum fluctuations, sometimes leading to 
highly non-classical behavior. Low-spin materials which 
constitute 'deformed' triangular antiferromagnets — that 
is, with spatially anisotropic exchange interactions that 
are not SU(2) symmetric — have indeed provided nu- 
merous surprises which are very likely quantum in ori- 
gin. Studies of the spin-1/2 compounds CS2CUCI4 and 
Cs2CuBr4 have been particularly fruitful in this regard. 
Inelastic neutron scattering has revealed striking dom- 
inance of a multi-par ticle c ontinuum in the dynamical 
response of Cs2CuCl j^^ l ^^ l This continuum is natu- 
rally explained in terms of spin-1/2 spinon excitations 
of weakly coupled chains. Ordered phases of this ma- 
terial show strong sensitivity to the magnitude and di- 
rection of the external magnetic fielcP^. The complex- 



ity of the h — T phase diagram has been attributed to 
the competition between several asymmetric exchanges of 
the Dzyaloshinskii-Moriya (DM) type and a weak inter- 
plane exchange interaction. In particular it has been 
suggestecP^ that, although weak, inter-plane exchange 
can dominate over stronger but frustrated inter-chain 
coupling and dictate the type of three-dimensional mag- 
netic ordering at low temperatures. 

Our study here is strongly motivated by the isostruc- 
tural material Cs2CuBr4, whose triangular planes are 
less anisotropic and exhibit weaker inter-plane cou- 
pling compared to CS2CUCI4, leading to rather differ- 
ent but equally rich phenomenology in magnetic fields. 
With fields directed in the plane of the triangular lay- 
ers, the experimentally determined phase diagram for 
Cs2CuBr4 hosts as many as nine phases^ at low tem- 
perature. A one-third magnetization plateau features 
very prominently'^HESl amongst the other, less under- 
stood, phases and offers a convenient starting point for 
theoretical analysia^S. Notably, CS2CUCI4 shows no signs 
of this plateau, which has been attributed^* to its more 
pronou nced s patial anisotropy and stronger inter-layer 
couplin^24|25j^ which is known to suppress the platearP^. 
And unlike the higher-spin material RbFe(Mo04)2 dis- 
cussed above^, the width of the magnetization plateau in 
Cs2CuBr4 is essentially T-independenl(2iJ, strongly hint- 
ing at its quantum origin.^ 23:2?] 

The sheer number of states present in the phase dia- 
gram of Cs2CuBr4, in comparison with the three phases 
expected from the standard isotropic Heisenberg antifer- 
romagnet on the triangular lattice, make it clear that a 
thorough theoretical study is required to begin under- 
standing this material's complex behavior. Here we at- 
tempt to access the global phase diagram of this sys- 
tem by analyzing the roles of three important known 
perturbations away from the 'ideal' classical triangular 
antiferromagnet — quantum fluctuations stemming from 
the low spin S = 1/2, spatial exchange anisotropy, and 
Dzyaloshinskii-Moriya (DM) interactions — all of which 
compete and favor different spin arrangements. Because 
the ideal classical model exhibits a large ground state 
degeneracy, we show that even weak spatial anisotropy 
and DM coupling are sufficient to qualitatively alter the 
standard Y-UUD-V phase diagram, stabilizing new spin 
orders including incommensurate 'umbrella' and planar 
states, a commensurate 'distorted V spin structure, and 
a commensurate 'inverted Y' phase. Further experiments 
(such as neutron scattering and NMR measurements) will 
be very helpful for identifying which of these orders ap- 
pear in the observed phase diagram, thereby sharpen- 
ing the outstanding theoretical questions that undoubt- 
edly remain. We also predict that — in sharp contrast to 
CS2CUCI4 — the phase diagram of Cs2CuBr4 ought to be 
insensitive to the direction of the magnetic field inside of 
the triangular planes. 

The remainder of the paper is organized as follows. We 
provide an overview of the model and strategy that we 
pursue in Sec. In) and review the physics of the isotropic 



triangular antiferromagnet in Sec. Ill Section IV ad- 
dresses the case where DM coupling is present but spa- 
tial exchange anisotropy is neglected. We then explore 
the Bose-Einstein condensation of magnons near the sat- 
uration field in Scc.lv) which allows us to simultaneously 
treat quantum fluctuations, DM coupling, and spatial 
anisotropy in a simple setting. The influence of spatial 
anisotropy on the global phase diagram is studied with 
and without DM interactions in Sees. VI and VII respec- 
tively. Finally, we provide a summary and concluding 
remarks in Sec. IVIIII 



II. MODEL AND STRATEGY 



We study the following Hamiltonian, 
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where the exchange integral Jr,r' is given by J on the 
horizontal bonds and J' on the diagonal zig-zag bonds as 
shown in Fig. )T) F or Cs2CuBr4, experiments have mea- 
sured the valuesP J = ll.iK and J' = 8.3K. The sec- 
ond term describes the Zeeman energy of spins in an ex- 
ternal magnetic field while the third, 7?dm, represents the 
asymmetric DM interaction between neighboring spins. 

Ideally, one would like to obtain the phase diagram 
for the quantum spin-1/2 problem above to begin un- 
derstanding the interesting phenomenology of Cs2CuBr4. 
In this paper we will attempt to access the physics of 
the spin-1/2 system coming from the large-S" limit, in- 
cluding quantum fluctuations perturbatively. This nev- 
ertheless still leaves a problem of substantial complex- 
ity, which can be understood by considering the classi- 
cal, spatially isotropic limit of H, without DM coupling. 
Let us denote this minimal classical Hamiltonian by Hq. 
As noted in the introduction and reviewed in detail in 
Sec. |III[ at finite magnetic fields Ho exhibits a large 'ac- 
cidental' ground state degeneracy. Quantum fluctuations 
(and thermal fluctuations at finite temperature) lift this 
ground state degeneracy, as do spatial anisotropy and 
DM coupling. However, all of these ingredients compete 
with one another, favoring completely different ordered 
states. To resolve this competition, our strategy will be to 
compute the lowest-order energy splittings for the degen- 
erate ground states of Hq coming from each effect. As we 
will see later quantum fluctuations and DM interactions 
split this degeneracy already at first order, while spa- 
tial anisotropy achieves this only at second order. Thus 
despite the fact that in Cs2CuBr4 spatial anisotropy [as 
quantified by {J — J') ~ 0.3 J] is expected to greatly ex- 
ceed the characteristic energy scales for DM interactions 
in the material, the two can in fact comparably influence 
the phase diagram. 

An additional complication arises from the fact that 
crystal symmetry of Cs2CuBr4 permits several DM 
term^i^ which together break the SU(2) spin symme- 
try enjoyed by the exchange coupling down to a discrete 



subgroup. One might then expect the phase diagram to 
depend sensitively both on the polar angle that h makes 
with respect to the z-axis, normal to the triangular plane, 
and the azimuthal angle h makes in the (x, y) plane. 
Such a highly anisotropic phase diagram indeed emerges 
in the more spatially anisotropic materiaP^ CS2CUCI4. 
In the limit of weak spatial anisotropy — which due to 
the relevance to Cs2CuBr4 is our main focus here — most 
of these DM couplings fortunately play an unimportant 
role. As justified in Appendix [Aj it indeed suffices to 
consider only 



ffDM = "5I^-[S"-X(Sr 
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where the DM vector is D = Dz (we assume D > 
throughout) and ^1^3 are vectors shown in Fig. flj The 
strength of this DM interaction is expected to be com- 
parable to that in the isostructural material CS2CUCI4, 
wherePLiR:! J/20. 

Notice that at h = the model with the above DM 
coupling still exhibits a U(l) spin symmetry correspond- 
ing to global spin rotations about the D vector. This 
leads to a major simplification — the phase diagram of 
H depends on the polar angle h makes with the z- 
direction but not on the azimuthal angle. All additional 
DM terms which lower this symmetry lift the acciden- 
tal degeneracy of Hq only at second order in their cou- 
plings, whereas Hum above has a first-order effect; see 
Appendix |X] for details. We emphasize that a nontriv- 
ial consistency check emerges here regarding relevance 
of our results to experiments on Cs2CuBr4. Our ap- 
proach postulates that the physics of this material can 
be accessed coming from the classical, isotropic model 
without DM coupling. If this postulate is correct, then 
contrary to CS2CUCI4, the experimental phase diagram of 
Cs2CuBr4 should be qualitatively insensitive to rotations 
of the field about the z direction. Recent experiments by 
Y. Takano and collaborators^^' have indeed shown that 
the low-temperature phase diagram of Cs2CuBr4 is qual- 
itatively the same for fields oriented along the material's 
b and c axes (x and y axes in our notation from Fig. 
Rj). This finding lends strong experimental support to 
our approach. 



III. REVIEW OF THE ISOTROPIC MODEL: 

J = J',D = 

A. Classical ground states 

In the absence of DM coupling, the classical isotropic 
Heisenberg model, where spins are described as 3D unit 
vectors, is well known to exhibit a large 'accidental' 
ground state degeneracy in applied magnetic fields*. The 
classical ground state structure can be conveniently ex- 
posed by expressing the Hamiltonian (up to a constant) 
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FIG. 1: Spatially anisotropic triangular lattice with ex- 
change J along horizontal bonds and J' along diagonal 
bonds. We define axes such that the sites lie in the (x, y) 
plane as shown. Vectors ^1,2,3 connect nearest-neighbor sites 
of the lattice. In the isotropic limit where J — J' and 
Dzyaloshinskii-Moriya coupling is absent, the Hamiltonian in 
Eq. ([l| exhibits a highly degenerate classical ground state 
manifold wherein spins order in an underconstrained three- 
sublattice pattern. We label the three sublattices by A,B, 
and C. 
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MA(r)=Sr + Sr+5i + S^+s, 
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represents the magnetization of an elementary triangle 
located at site r. Below the saturation field /ig^t — 9 J, it 
is clear from Eq. pi that classical ground states satisfy 
the condition MA(r) — h/(3J) for all r. States fulfilling 
this requirement exhibit a three-sublattice structure as 
shown in Fig. Ill whose spins Sa.b,c are determined from 



Sr7 = 



3J' 



(5) 



The ground state degeneracy follows immediately from 
Eq. ([5]) , since the spins on the three sublattices are speci- 
fied by a total of six angles which are constrained by only 
three equations. Parametrizing the spins on sublattice a 
as 

Sq = (sin 9a cos (pa , sin Oa sin (/>„ , cos 0a), (6) 

the angles must specifically satisfy 

h • X 

sin 9ji cos (f>A + sin Og cos cpg + sin 9c cos (jjc 

sin 9j^ sin 0^ + sin 9b sin (jjg + sin 9c sin c/jc 
cos 9 a + cos 9 b + cos 9c 



3J 
h-y 

3J 
h z 

I3X 



(7) 



Equations ([7]) allow one to, say, express (f>B,4'Cj and 9c 
in terms of (/)a, 9a, and 9b- Note, however, that classical 
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FIG. 2: Important members of the classical ground 

state manifold exhibited by Eq. ([T| when J = J' and 
Dzyaloshinskii-Moriya coupling is absent: (a) Y, (b) up-up- 
down (UUD), (c) V, (d) 'distorted V, (e) 'inverted Y', and 
(f) umbrella states. The arrows denote the spin orientations 
on the A, B, and C sublattices of Fig. fT] 



ground states do not exist for all possible values of the 
latter angles. 

At h = where the classical Hamiltonian exhibits 
0(3) symmetry, this degeneracy is symmetry-related: the 
three unconstrained angles reflect the arbitrariness of the 
plane in which the 120° spiral orients and the freedom 
for rotating all spins by an arbitrary angle within that 
plane. Introducing a finite magnetic field reduces the spin 
symmetry down to U(l). The three free angles neverthe- 
less persist, only one of which is now symmetry- related 
(corresponding to global rotations of the spins about the 
magnetic field axis). Consequently, in finite magnetic 
fields below saturation, Hq exhibits an accidental classi- 
cal ground state degeneracy specified by two continuously 
deformable angles. Figure |2] illustrates some important 
members of this ground state manifold which we will fre- 
quently refer to later on: (a) Y, (b) up- up-down (UUD), 
(c) V, (d) 'distorted V, (e) 'inverted Y' , and (f) umbrella 
states. 



B. Symmetries of the classical ground states 

It is important to emphasize the symmetry distinction 
between the above phases. Because (Sr • h) is uniform in 
umbrella configurations, this order is exceptional among 
the classical ground states in that it preserves the original 
lattice translation symmetry [when followed by an ap- 
propriate global U(l) spin rotation]. All other classical 
ground states exhibit a nontrivial three-sublattice pat- 
tern of (Sr'h) and therefore spontaneously break discrete 
translation symmetry. UUD order is also exceptional — it 
alone preserves global U(l) spin rotations about the field 
axis and therefore does not break any continuous symme- 
try. Ultimately, the U(l) spin symmetry enjoyed by the 
UUD state permits the formation of a one-third magneti- 



zation plateau in this system. In a boson representation 
of the spins, umbrella order corresponds to a superfluid 
phase, UUD order is a solid phase, and all other classical 
ground states constitute supersolids.'SSEil 

These 'supersolids' can be further distinguished by 
symmetry. Since in the V state spins on two sublattices, 
say A and B, are parallel, this order is symmetric under 
spatial 7r/3 rotations about a site on the C sublattice. 
This can also be phrased as invariance of the V state 
with respect to permutation of the A and B sublattices. 
The distorted V state, while smoothly connected to the 
V state, violates this symmetry (the A and B sublattices 
are not equivalent anymore) and therefore constitutes a 
distinguishable phase in the isotropic system. Finally, in 
both the Y and inverted-Y states spins on one sublattice, 
say C, point either along or against the field, and both 
orders are invariant under 7r/3 rotations about a site on 
sublattice C followed by a global tt spin rotation. While 
the latter phases exhibit identical symmetries, one can 
not deform the spins smoothly from one configuration to 
the other without breaking additional symmetries; thus 
it is still meaningful to classify them as different states. 

All of the orders displayed in Fig. [2] thus consti- 
tute distinct phases when the full symmetries of the 
isotropic Heisenberg triangular antiferromagnet in field 
are present. Note, however, that when the symme- 
try is lowered by including DM coupling and/or spatial 
anisotropy, this ceases to be the case. The V and dis- 
torted V states, for example, are not distinguishable once 
7r/3 rotation symmetry is removed by adding either of 
these ingredients. Similarly, the distinction between the 
UUD and Y states hinges on U(l) spin symmetry and 
rotation symmetry, both of which are lost in the pres- 
ence of DM coupling when the magnetic field has a finite 
component in the (x, y) plane. The UUD and distorted V 
states can, by contrast, remain distinct even in this phys- 
ical situation, despite the very low symmetry remaining 
for the problem. For example, when h = /ix and DM 
coupling is non-zero the Hamiltonian remains invariant 
under S^ -^ S% and S'^>^ -^ -Sl'^; the UUD state re- 
spects this discrete symmetry while the distorted V state 
does not. We will return to these issues in subsequent 
sections of the paper. 



C. Influence of thermal fluctuations 

Thermal fluctuations provide one mechanism that lifts 
the large accidental ground state degeneracy of the classi- 
cal model. At finite temperature the system minimizes its 
free energy, F — U — TS, and although the ground states 
exhibit identical energies U, their entropies generically 
differ. Thus only the most entropically favored states 
emerge. The finite temperature phase diagram was es- 
tablished long agc)^, and is reproduced with our numerics 
in Fig. [3} 

This set of Monte-Carlo simulations is based on ALPS 
code.l^ System sizes ranged from 24 x 24 to 96 x 96 
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FIG. 3: (Color online) Phase diagram of the classical isotropic 
model. Amongst the large set of accidentally degenerate clas- 
sical ground states, thermal fluctuations select the Y, UUD, 
and V states illustrated in Figs. pTa) through (c). 



with most simulations being carried out on 48 x 48 lat- 
tices. (The triangular lattice can be viewed as a square 
lattice with bonds along one diagonal; the system sizes 
above refer to the dimensions of such a square lattice. 
All simulations were performed with periodic boundary 
conditions.) Every (T, h) coordinate was simulated using 
the standard Metropolis algorithm, with 500,000 Monte 
Carlo steps for thermalization and another 500,000 for 
measurement in the 24 x 24 system. For 48 x 48 lattices 
these numbers were in the range 150, 000 — 250, 000. All 
Monte Carlo data presented here and below correspond 
to J = 1. Different phases and boundaries between them 
have been identified by the behavior of the specific heat, 
the magnetization M and its field derivative dM/dh, and 
the vector chirality k and coplanarity K defined below. 
The vector chirality is defined as 
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where N is the total number of spins in the system and 
the normalization factor ensures that the maximal chiral- 
ity magnitude is unity. Taking h. — hz for concreteness in 
the remainder of this section, it is useful to consider the 
longitudinal k^ and transverse k± — (k^-|-k^)^'^ compo- 
nents of the vector chirality, in addition to its magnitude 
K = {kI + K^y^'^. Finite longitudinal chirality k^ identi- 
fies non-coplanar spin structures of umbrella type while 
a non-zero transverse chirality signifies coplanar ordering 
in which no two spins in the unit cell are parallel (such 
as Y states). Notice that by construction the chirality 
vanishes in three-sublattice states containing two paral- 
lel spins, such as the UUD and V states. 

The coplanarity K, introduced in Ref. 1331 provides 



another useful indicator of planar spin structures. This 
quantity follows from the A/B/C sublattice magnetiza- 
tions 
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by constructing 



Kab = (Ma X Mb) x h 



(9) 
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and similarly for K-bc ^^nd H-cA ■ The coplanarity is then 
given by the combination 
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As defined K is finite for any coplanar spin configuration 
which includes the magnetic field in its plane. In particu- 
lar, K is finite in the V state but vanishes in the collinear 
UUD configuration. 

Our Monte-Carlo simulat ions , summarized in Fig. |3J 
agree well with existing datcP^ and clearly demonstrate 
entropic selection of the Y, UUD, and V states over non- 
coplanar order. This is consistent with the usual intuition 
that entropy disfavors non-coplanar structures a nd fa- 
vors coplanar and, when available, collinear orderingl^Iil]. 
Coplanarity of the selected states is reflected in an ex- 
ceedingly small value of the longitudinal chirality Kz over 
the e ntir e ma gnetic field range < h < 9J. Figures 
4(a) and |4(b)] exhibit the magnetization A/ and its field 
derivative, dM/dh, a.tT — 0.2 J over the same field range. 
Notice that the UUD state is very well identified by the 
abrupt variation of dM/dh, with two peaks in Fig. |4(b)| 
representing the lower and upper critical fields. Away 
from the UUD state the slope of the magnetization ap- 
proaches 1/(9J), which is just the uniform susceptibil- 
ity of the classical antiferromagnet at zero temperature. 
The magnetization slope in the UUD state [Fig. 4(a)| 
is visibly smaller, but clearly remains finite. This is a 
manifestation of the 'dual' role played by temperature: 
it simultaneously selects coplanar states and thermally 
disorders them, eventually resulting in a paramagnetic 
state above a field-dependent critical temperature. 

Despite being such a well-studied problem, there re- 
mains a serious question about the phase diagram of the 
isotropic triangular lattice antiferromagnet at h = 0. As 
described above in Sec. [Ill A[ the large degeneracy of 
the model in this limit is symmetry related and refiects 
arbitrariness of the ordering plane in which 120° spiral 
forms. This makes the order parameter space isomorphic 
to SO (3), the group of rotations of a three-dimensional 
rigid bodySl. According to the Mermin- Wagner theorem, 
this continuous symmetry cannot be broken at any finite 
temperature. Nonetheless Monte-Carlo simulations, in- 
cluding ours, do show a weak peak in specific heat at 
finite ternperature, approximately 0.33J. It has been 
suggestecPthat this finite-T feature in fact reflects a non- 
trivial phase transition associated with binding of Z2 vor- 
tices permitted by the S0(3) structure. At present our 



study, which aims to understand the global features of the 
deformed Heisenberg model at finite T and h, has noth- 
ing to add to this interesting issue^^HMl^ ^j^g resolution of 
which requires more extensive Monte-Carlo simulations. 

We should emphasize, however, that at finite mag- 
netic fields the entropically selected states break trans- 
lational symmetry and introduce three inequivalent sub- 
lattices A, B and C. Consequently, at any finite h up 
to saturation there exists a finite critical temperature Td 
arising from discrete symmetry breaking associated with 
the spontaneous selection of the sublattice ordering. In 
Fig. [3] this hull-shaped Td{h) separates low-temperature 
ordered states from a high-temperature paramagnetic 
phase in which only the uniform magnetization M is 
present. 

The finite-field situation contains yet another unre- 
solved issue, related to the U(l) symmetry correspond- 
ing to rotations about the magnetic field direction. The 
coplanar Y and V states break this symmetry by se- 
lecting an ordering direction in the plane perpendicular 
to the field direction. While due to dimensionality the 
U(l) symmetry can not be truly broken at finite tem- 
perature, there exists a Kosterlitz-Thouless (KT) tem- 
perature Tkt below which the transverse magnetic order 
exhibits power-law correlations. In principle one could 
have Tkt 7^ T^ (see Ref. 37 for a recent example of such 
behavior). From the available data, however, it appears 
that within our numerical accuracy the two transitions 
coincide, Tkt = 2d. A related issue is that at finite tem- 
perature, the transitions from UUD to either Y or V are 
also of KT type, and as such are inherently broad and 
difficult to locate precisely in numerica^Sl. We would like 
to emphasize again that our main focus is on identifying 
possible phases rather than studying the nature of the 
transitions between them. Readers interested in prop- 
erties of such phase transitions in a similar context are 
instead referred to the recent detailed study in Ref. l40l 
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FIG. 4: (Color online) (a) Magnetization M and (b) its field 
derivative dM/dh in the vicinity of one-third magnetization. 
The data were obtained with Monte Carlo simulations of 
the classical isotropic model at temperature T = 0.2J. The 
coUinear UUD state underlies the drop in magnetic suscepti- 
bility near h — 2.5J, visible in both (a) and (b). The dashed 
line in (a) shows the T — magnetization of the classical 
triangular antiferromagnet, M = h/{9J). 



As first shown by Chubukov and Goloso'v^, the acciden- 
tal degeneracy inherent in the classical isotropic model 
can be lifted even at zero temperature by quantum fluc- 
tuations. Here zero-point energy of spin waves — rather 
than entropy — provides the degeneracy-lifting mecha- 
nism. Quantum ground state selection can be system- 
atically explored using the machinery of the 1/5 expan- 
sion. To this end, one starts with a particular clas- 
sical ground state in the degenerate manifold and de- 
fines rotated spin coordinates 8^/5/0(1") on the three 
sublattices such that the classical order corresponds to 



^A/B/C 



(r)) 



Si. One then introduces Holstein- 



Primakoff bosons a(r),&(r), and c(r) to express the ro- 



tated spin operators to leading order in 1/S as 



Slir) = y|[a(r)-|-at(r)] 






n/-[a(r) 



at(r)] 



= S-a\r)a{r) 



(12) 



and similarly for Ss/c(i")- Retaining only the leading- 
order terms in 1/S generates a Hamiltonian of the form 
H = Eq + H2, where Eq is the classical ground state 
energy and H2 is quadratic in the three boson fields. 
Although at a given magnetic field /i, Eq is identical for 
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FIG. 5: (Color online) Leading 1/ S quantum energy splittings 
between the Y/UUD/V states and the inverted-Y (dashed 
line) and umbrella states (solid line). These curves illustrate 
an important trend — quantum fluctuations are most effective 
at splitting the accidental degeneracy near 1/3 magnetization 
where low-energy coUinear states are available. 



all allowed classical ground states satisfying Eqs. (l7|, the 
spin-wave Hamiltonian H2 does depend on the particular 
ordered configuration around which one is expanding, via 
the dependence of the spin- wave dispersion lo^ on the an- 
gles Qajbjc ^-iid 4>A/B/c- The zero-point energy for spin- 
waves, {H2) = (1/2) X];.a;fc, therefore yields a 1/5* cor- 
rection to the energy that lifts the classical degeneracy. 
See, for example, Ref. 41 for an earlier application of this 
strategy to a related problem. We note that this analy- 
sis can be carried out rather efficiently by parametrizing 
the original spin operators via Eq. ([6]) and deriving the 
general form of H2 as outlined above. The constraints in 
Eqs. ([7]) allow one to express (f)B^ 4'c, and 6c in terms 
of (j)A, Oa, and 9b] due to the U(l) spin symmetry that 
is present when h 7^ 0, one can arbitrarily set (J)a — 0. 
Consequently, H2 can be expressed in terms of only two 
angles, 6a/b- Diagonalizing H2 for a mesh of 9a and 
9b then allows one to compute the zero-point energies 
for the classical ground state manifold to deduce which 
states quantum fluctuations favor. Although the selec- 
tion for the isotropic model is already well-known, we 
will use this scheme later on to explore the competition 
between quantum fluctuations and DM interactions. 

It is important to emphasize a subtle point in the pro- 
cedure outlined above. At order 1/S quantum fluctua- 
tions not only split the degeneracy through zero-point 
motion, but also renormalize the spin structure in the 
classical ground states. Thus to order 1/S the spin an- 
gles should read 9a = 9^ + 69a , and similarly for other 
angles, where 9^ corresponds to a classical ground state 
and 69 A is a 1/S quantum correction (which we do not 
calculate here). When evaluating the zero-point ener- 
gies (7^2) to order S, clearly one can safely neglect such 
quantum corrections to the spin states. One may worry. 



however, that since the classical ground state energy Eq 
is proportional to 5*^, evaluating Eq with the renormal- 
ized spin angles may produce an order-5 contribution to 
the energy that is comparable to the leading zero-point 
energy. This is certainly not the case — because the clas- 
sical ground states minimize Eq, the shifts 69a appear 
only at second order and therefore contribute to the en- 
ergy only at order 5"; see, e.g., the discussion in Ref. 
Hni Thus while renormalization of the spin structure 
is crucial for capturing certain physical quantities such 
as quantum corrections to the magnetization curve as a 
function of field, this effect is indeed unimportant for our 
main purpose: evaluating the ground state energies to 
leading order in 1/S. 

As in the case of thermal fluctuations, quantum fluc- 
tuations disfavor non-coplanarity and prefer collinear 
order, leading to Y states below 1/3 magnetization 
{0 < h < /isat/3), UUD ordering at 1/3 magnetiza- 
tion {h = ft.sat/3), and V states at larger magnetiza- 
tion (/isat/3 < h < ft-sat) up to saturatiorP. While it 
is well known that the UUD phase realizes a magneti- 
zation plateau also in the quantum problerrP^ we em- 
phasize that this fact arises from higher-order quantum 
corrections such as the 5"° effect discussed above. One 
can recover this feature by incorporating 1/S quantum 
renormalization of the equilibrium angles for the Y and 
V states, which broadens the UUD state into a plateatP. 
This result can be also be seen in the following comple- 
mentary way. Since the UUD state is collinear it does not 
break any continuous symmetry, and therefore exhibits 
only gapped spin excitations when one renormalizes the 
spin-wave spectrum to leading order in 1/5' (note that 
such a renormalization shifts the ground state energy 
only at order S^). Thus a finite change in the magnetic 
field is required to destabilize this phase. 

To establish intuition that will prove beneficial in sub- 
sequent sections, it is useful to now ask how robust the 
quantum selection of Y/UUD/V order is at different mag- 
netic fields. One can quantify this by computing the first 
1/S quantum correction to the energy, i.e., evaluating 
{H2), for the Y/UUD/V states and comparing it with the 
corresponding correction for other classical ground states 
which quantum fluctuations disfavor. The size of these 
energy differences AE = (1/2) ^Ja;^*'^" - uj^^'"™^^] 
provide one with a rough guide for how robust the quan- 
tum ground state selection is to the inclusion of other per- 
turbations, such as DM interactions. Figure [5] presents 
the energy splittings AE obtained for the umbrella (solid 
line) and inverted-Y states (dashed line) as a function 
of h. To produce the data shown we have extrapolated 
large-S* corrections to 5* = 1/2. It is seen that Y/UUD/V 
order is preferred over umbrella and inverted-Y states for 
all values of the field. It is also seen that being co-planar, 
the inverted-Y state has lower energy than the umbrella 
one. Although only two states, umbrella and inverted- 
Y, are considered, the figure clearly illustrates a gen- 
eral trend which recurs throughout the paper — quantum 
fluctuations are most effective at selecting ground states 



in the vicinity of 1/3 magnetization where low- energy 
collinear spin configurations are available. 



E. Biquadratic approximation of quantum 
fluctuations (zero-point motion) 

There is an alternative, albeit less rigorous, way to ac- 
count for the influence of quantum fluctuations, within 
a purely classical spin Hamiltonian. This approach 
is based on the observation that coplanar and espe- 
cially collinear spin states — which quantum fluctuations 
tend to favor most strongly — can be stabilized classically 
upon addinga suitable biquadratic spin coupling to the 
HamiltoniarP. This biquadratic interaction represents 
an effective Hamiltonian which is generated by quantum 
and/or thermal fluctuations"*"^ . We have incorporated 
such a biquadratic spin interaction and optimized its 
(field dependent) coupling constant such that we repro- 
duce the leading 1/S energy differences between several 
important states in the classical ground state manifold. 
In this way we find that quantum zero-point energy can 
be semi-quantitatively described by the following classi- 
cal biquadratic Hamiltonian with negative coupling con- 
stant. 



5Hu = -A{h)J2iSr-Sr')\ 

(rr') 

A/,N 0.0268 
Aih) = -^ 



(13) 



f - 0.03/1 V/^sat - h 



Note that spins Sr appearing in dH\^i are classical unit 
vectors and that all dependence on the magnitude of the 
microscopic spin S is contained in the field-dependent 
coupling constant A(/i). Figure p] compares the en- 
ergy differences between various classical ground states 
specified in the caption using the leading 1/S results 
(solid lines) and the biquadratic classical approxima- 
tion (dashed lines), showing the quite good quantitative 
agreement between the two approaches. 

The full benefit of approximating quantum zero- 
point energy by the classical biquadratic interaction will 
become clear when we address the effects of spatial 



anisotropy J' 7^ J in Sees. VI B and VHB The clas- 
sical energy gain due to spatial anisotropy is quadratic 
in the difference J—J'. Such second-order corrections are 
not in general easy to account for analytically, making it 
difficult to consistently compare the anisotropy-induced 
energy splittings with those arising from quantum ef- 
fects and DM interactions. Using the biquadratic inter- 
action in Eq. (131, however, allows us to circumvent this 



problem and find the result of the competition between 
the different perturbations by numerically studying, via 
standard Monte Carlo techniques, the classical system 
described by Hq + i?DM + SHhi- Later we will refer to 
ground states obtained within such a classical model with 
biquadratic spin couplings as 'pseudo-quantum ground 
states' to emphasize that despite corresponding to classi- 
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FIG. 6: (Color online) Energy differences between vari- 
ous states using the leading 1/5 result from spin-wave the- 
ory (solid lines) and the biquadratic classical approximation 
(dashed lines). Green lines represent the difference between 
Y/UUD/V and inverted Y energies; red between Y/UUD/V 
and umbrella; and black between inverted Y and umbrella. 



cal spins, such states are expected to semi-quantitatively 
reflect the influence of quantum fluctuations. 

As a warm-up exercise of such a strategy, we first 
simulate a simpler "isotropic plus biquadratic" Hamil- 
tonian iJiso+bi = Ho + SHbi. Figure [t] shows dM/dh 
and the coplanarity K for this model at a low tempera- 
ture of T = 0.027 J. This low temperature was reached 
via simulated annealing to avoid becoming trapped in 
local (but not global) free-energy minima. Specifically, 
Monte Carlo simulations were performed on a 24 x 24 
system starting at T = 0.4J, then cooling down in small 
AT = 0.009J steps. At each temperature the system 
was equilibrated with 35,000 Monte Carlo steps. Mea- 
surements were taken during the last 30,000 steps. 

The pronounced minimum of dM/dh visible in Fig. [t] 
over the field interval 2.2J < /i < 3.2J identifies this 
region with the collinear UUD state. The coplanar na- 
ture of the adjacent phases is evident from the plot of 
K vs. h in the same Figure. Note that despite the very 
low temperature considered, a substantial UUD plateau 
remains whose width greatly exceeds that of the purely 
classical model without interactions; see Fig. [3] for com- 
parison. The width of the plateau is in fact proportional 
to A(/i) as discussed in Ref. I9l quantum fluctuations, 
modeled here by SH^i, produce a stable UUD magnetiza- 
tion plateau over a finite field interval, even at T = 0. We 
note that the usual spin- wave approach predicts a plateau 
over a field interval of width Ah « 1.8J/{2S), which is 
somewhat larger than the plateau width obtained in our 
isotropic-plus-biquadratic simulations. 
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FIG. 7: (Color online) dM/dh (blue circles) and coplanarity 
K (green squares), defined in Eq. (Ill, for the classical 
"isotropic plus biquadratic" model at T = 0.027 J, obtained 
by annealing from high temperature. 



IV. SPATIALLY ISOTROPIC MODEL WITH 
DZYALOSHINSKII-MORIYA INTERACTIONS: 

J = J',D/0 

In this section we will explore how DM coupling af- 
fects the classical and quantum phase diagrams discussed 
above. One complication that arises here is that once 
we invoke DM interactions, the magnetic field direction 
is no longer arbitrary as it was in the previous section. 
Recall that the DM vector points along the z-axis. We 
will begin in Sec. |IVA| with the case where the field is 
applied along the z-direction, perpendicular to the tri- 
angular lattice plane. This field orientation is simplest 
to analyze since the Hamiltonian including Hbm then 
preserves global U(l) spin symmetry about the z-axis. 



In-plane field orientations are studied in Sec. IV B 



We should also emphasize that in any field orientation, 
DM coupling breaks 7r/3 rotation symmetry since i?DM 
in Eq. (pi) only couples spins on neighboring 'chains' of 
the lattice. Thus the model with only DM coupling (and 
no spatial anisotropy) is clearly fine-tuned. Neverthe- 
less, by examining the effects of DM alone we will gain 
valuable intuition that will be useful when we include the 
additional complication of spatial anisotropy in Sec. |VII 



A. Parallel orientation: D || h = /iz 

1. Classical Ground States 

We will first explore how DM interactions lift the ac- 
cidental degeneracy of the classical isotropic model using 
first-order perturbation theory. At this order, one simply 
needs to evaluate the DM energy using the unperturbed 
form of the classical ground states found in the previous 
section. Using the parametrization for the ground states 



in Eq. (l6| and the constraints in Eqs. (I?]), we obtain 



(iJoM) = 2(jND{ - sin^ 6Ia sin^ 6Ib 



1 



324J4 



[h^ + 9J^ - 6Jh{cos9A + cos0b) 



18 J^ cos 6'a cos 6'b] 



1/2 



(14) 



where a = sign [sin (^^ — 4'b)] determines the chirality 
of the spins in the (x, y) plane. The DM energy above 
is minimized when cos 6*^ — cos 9b — ^ and cr = — 1. 
One can readily see using Eqs. ([7]) that this corresponds 
to umbrella states of Fig. I2jf ) for all fields (with a spe- 
cific chirality because D ^ breaks inversion symme- 
try). This outcome is quite natural given that DM cou- 
pling at h = favors coplanar spiral order with spins 
pointing in the triangular lattice plane; the field simply 
cants the spiral out of the plane, producing an umbrella 
pattern. Note that this is markedly different from the 
coplanar and coUinear configurations favored by thermal 
and quantum fluctuations. This leads to an interesting 
competition once one incorporates these additional ingre- 
dients, as we now address. 



2. Influence of thermal fluctuations 

One can deduce the qualitative outcome of the compe- 
tition between DM interactions and thermal fluctuations 
at flnite temperature using the following general argu- 
ment. In the previous section we discussed that while at 
T — the isotropic model exhibits many classically de- 
generate ground states, at flnite temperature entropy lifts 
this degeneracy in favor of Y/UUD/V states. Moreover, 
this selection is strongest at intermediate temperatures 
where thermal fluctuations provide a large free-energy 
splitting but are not so severe as to destroy ordering alto- 
gether. (This can be seen most clearly through the width 
of the UUD plateau in Fig. [3J which reaches at maximum 
at around T = 0.35J before entering the paramagnetic 
phase at slightly higher T.) 

Now suppose one slowly turns on DM interactions. At 
zero temperature any flnite DM coupling will immedi- 
ately select umbrella ordering since the entropic split- 
tings vanish as T — ^ 0. By continuity umbrella states 
must persist over a flnite temperature window that de- 
pends on both fleld and the DM coupling strength. Sim- 
ilarly, at intermediate temperatures the entropically fa- 
vored Y/UUD/V states must also by continuity survive 
some amount of DM coupling. This leads to multi-stage 
ordering as a function of temperature. Consider for illus- 
tration the interesting example of one-third magnetiza- 
tion with weak DM coupling. As temperature increases 
from zero the system will begin in an umbrella state, 
transition into an UUD plateau at intermediate temper- 
atures, and eventually give way to a paramagnet at high 
temperatures. 
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energies for a dense subset of classical ground states as 
described in Sec. 



HID 



one can determine the spin order 
selected when both competing effects are present. 

Carrying out this procedure with D/J = 0.05 and 
S = 1/2, we have found that the umbreUa phase is the 
quantum ground state for sufficiently small magnetiza- 
tion {h/hss^t ^ 0.3) and, also, sufficiently near saturation 
{h/hsa.t ^ 0.94); see Fig. [sja) which summarizes our nu- 
merics for this case. In between, however, a different 
ground state emerges which reflects a compromise be- 
tween DM interactions and quantum fluctuations — the 
distorted V phase of Fig. ^d). This state resembles the 

Y and V orders favored by quantum fluctuations, but is 
distorted in a non-coplanar fashion to gain DM energy. 
Note that the distortion of the spins away from perfect 

Y or V order is not unexpected: these states, as well 
as UUD, exhibit rotation symmetry which is explicitly 
broken by the DM interaction in Eq. (pi) (recall the dis- 
cussion at the end of Sec. IlIB ). The analytical structure 
of the distorted V state is easy to obtain in the high 
magnetic field limit as described in Sec. |VB[ 



FIG. 8: (Color online) Phase diagram of the quantum Heisen- 
berg model with DM interactions, as a function of magnetic 
field h. The magnetic field orients along the DM axis, D || h, 
in (a) and orthogonal to the DM axis, D _L h, in (b). 



B. Orthogonal orientation: D _L h = /ix 



Classical Ground States 



The qualitative picture that emerges is that as one 
increases DM coupling from zero, DM interactions eat 
away at the entropically-stabilized regions from the low- 
temperature side, eventually removing the Y/UUD/V 
phases entirely in favor of umbrella order at sufficiently 
large coupling. Obtaining a more quantitative under- 
standing of this competition requires extensive Monte 
Carlo simulations, which we will not explore here. We 
will instead perform such a study in the more physi- 
cally relevant case with spatial anisotropy in Sec. |VI A[ 
where very similar physics arises. It should also be kept 
in mind that in this subsection we neglected quantum 
effects entirely. Below we incorporate quantum fluctua- 
tions and show that they lead to a substantially richer 
zero-temperature phase diagram than that of the purely 
classical model. 



3. Quantum ground states 

To incorporate DM interactions and quantum fluctu- 
ations at T = 0, we treat 1/5 and D/J as expansion 
parameters of the same order of magnitude. The first- 
order corrections to the classical ground state energies 
can then be obtained simply by adding the 1/S spin- wave 
contribution discussed in Sec. IIIIDI to the classical DM 
energy in Eq. (JMJ) . [For example, the effect of quantum 
fluctuations on the DM energy produces a higher-order 
correction ^ {D/J){l/S).] By evaluating the perturbed 



We will now extend the analysis of the previous sec- 
tion to the case where the magnetic held is applied in the 
plane of the triangular layers, along the x-direction for 
concreteness. This field orientation is complicated by the 
fact that the quantum Hamiltonian including DM cou- 
pling now lacks U(l) spin symmetry; instead the Hamil- 
tonian is symmetric only under discrete spin rotations 
S^ -^ S% and S*^^^ -^ -3^4. As above, we start by 
deducing the order favored by DM coupling in this field 
orientation using first-order perturbation theory. Using 
Eq. ^ to parametrize the classical ground states along 
with the constraints in Eqs. Q, the first-order correction 
to the energy arising from DM coupling can be written 

(-^dm) = — — [/i(sin6lBsin0B - sin^Asin^^) 

-|- 9Jsin6'^sin0B sin((?!)^ — (/)b)], (15) 

where the angles are subject to the additional constraint 
h 



1 



h 
3J 



3J 



(sin 9a cos (f>A + sin Ob cos (Jjb) 



sm a A sm ub cosyipA 
cos 9a cos 9b ■ 



(16) 



One can use Eq. ( 16 ) to eliminate 9b from Eq. ( 15 ) and 



then minimize the DM energy over the remaining angles 
9A^4'A^ and (pB- In this way we find that DM coupling 
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selects coplanar inverted-Y states of Fig.[2je), with 

Oa/b/c = ""/S 



(pA = -(PB = cos 

<^c = 0. 



1 

2 V3J 



- 1 



(17) 



The emergence of planar ground states for this field 
orientation is quite natural, again because DM interac- 
tions favor spins orienting within the triangular lattice 
plane. Among the various planar arrangements, inverted- 
Y states most effectively gain DM energy by keeping 
neighboring spins as far from parallel as possible (within 
the constraints set by the ground state manifold) . In con- 
trast, the UUD and V states favored by quantum fluc- 
tuations yield a vanishing first-order DM energy because 
in both configurations spins on two sublattices are par- 
allel. We note that the favorability of inverted-Y over Y 
states at low fields is less obvious; the latter also gains 
DM energy, though to a lesser extent than the former. 

Just as in the h. = hi field orientation, incorporat- 
ing thermal or quantum fluctuations therefore leads to a 
delicate competition with DM interactions. The discus- 



sion from Sec. IV A 2 regarding entropic effects applies in 
this field orientation as well. Perhaps the most impor- 
tant qualitative distinction with the h. ^ hi field orien- 
tation is the additional loss of U(l) symmetry here due 
to the orthogonality between h and D vectors. This im- 
plies the demise of the 1/3 magnetization plateau since 
the UUD and Y states then no longer constitute distinct 
phases separated by a spontaneous U(l) symmetry break- 
ing. Sec. |VIIB] discusses this effect in greater detail; see 
Fig. [15] for an illustration of the rounding-off of the mag- 
netization plateau by DM interactions. 

In the following subsection we will resolve the outcome 
of the competition between quantum fluctuations and 
DM interactions at zero temperature, which is qualita- 
tively different from the previously considered h. = hi 
case. 



2. Quantum Ground States 

The quantum zero-temperature phase diagram in this 
field orientation is obtained as described in Sec. |IV A 3[ 
again treating quantum corrections and DM coupling 
to first order with D/J — 0.05 and S extrapolated to 
1/2. Figure [8](b) depicts our results. Quantum effects 
and DM coupling compete in a subtle manner for this 
field orientation. The former dominates at low fields, 
where Y states prevail, while the latter dominates at 
high fields where inverted-Y states appear. Distorted 
V order, which reflects a nontrivial compromise between 
the two competing effects, appears at intermediate mag- 
netization, 1/3 < h/hsat ^ 0.65. Note that all states 
are either perfectly planar or very nearly so here, which 
again is rather natural. Also notable is the absence of 
the collinear UUD state in the quantum phase diagram. 



As argued above, this is caused by the loss of U(l) spin 
symmetry in the D _L h geometry. 

Observe that the quantum phase diagrams predicted 
for the two field orientations, h oc z and h ex x, are 
very different. At either low or high fields, this implies 
that at least one phase transition separating the umbrella 
and Y/inverted-Y states ought to appear as one rotates 
the field from the z-direction to the x-direction. Recent 
experiments^ do indeed find a very complicated evolu- 
tion of the phase diagram as the field rotates in this plane, 
featuring several intervening phase transitions. General- 
izing the results of this paper to explore this crossover 
would be an interesting future research direction. 



V. BEC ANALYSIS NEAR SATURATION 

We will now shift our focus to phases realized at high 
magnetic fields near saturation, which can be elegantly 
described via Bose-Einstein condensation of magnon ex- 
citations above the fully polarized (saturated) state^SHU, 
In this way magnetically ordered phases just below sat- 
uration {h = hsut — Sh) can be understood as a BEC 
instability of the ground state just above saturation 
{h — hsat + 5h). This limit has the advantage of al- 
lowing one to incorporate all of the competing effects 
of interest — quantum fluctuations, DM coupling, and 
spatial anisotropy — in a single controlled and coherent 
framework. 



A. Preliminaries: J = J' ,D — 

At h > hsat the ground state of the Heisenberg Hamil- 
tonian is known exactly: it is given by the fully polarized 
eigenstate of the Hamiltonian in which all spins point 
'up'. The excitations are magnons, i.e., plane- wave states 
of overturned spins. These are conveniently described by 
Holstein-Primakoff bosons, similar to Eqs. (12). Since all 



spins point in the same direction, it suffices to consider 
only a single boson species here; thus we write 



S^{r) = S- a^(r)a(r), S+{r) = V2S a(r), (18) 

for all r. The isotropic Heisenberg Hamiltonian with- 
out DM coupling then reads, neglecting 1/S and smaller 
contributions, 



H 

Ho 



Ho + V 

Y, 5[J(k) - J(Q) - f,]ala^ 



^ = i^ 5][2J(q)-J(k)-./(k'-q)] 



k,k',q 



k+q"^k'-q 



ak'flk- 



(19) 
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Here 

3 

J(k) = y^ Jr,r+a, cos(k- Sj) 

= J [cos(fc^) + 2 cos{kx/2) cos{V3ky/2) (20) 
is the Fourier transform of the exchange interaction, and 



component is modulated with wavevector 2Q. The copla- 



M = (^sat - h)/S 



(21) 



is the chemical potential which controls the state of the 
bosonic system. For h > /igat the chemical potential is 
negative and the ground state is the boson vacuum |0) 
satisfying ak|0) = for all k. The wavevector Q in 
Eqs. ( 19 ) corresponds to the minimum of the magnon 

mink[JYk)]. A pecu- 
and indeed 



energy and follows from J(Q) 

liarity of the triangular antiferromagnet^, 

many other frustrated spin system^^, is the fact that 

there are two distinct wavevectors minimizing J(k). For 

the isotropic two-dimensional triangular antiferromagnet 

these are ±Q — (±47r/3, 0). Hence, at h = hgat (or 

equivalently fi — 0), the magnon dispersion touches zero 

simultaneously at ±Q. One then needs to understand 

whether it is energetically favorable to condense magnons 

at one or both of these wavevectors, i.e., form a single-Q 

or double-Q condensate. 

The analysis proceeds by parametrizing the magnon 
modes as follows. 



Ok 



Ntp+qSk^Q + vN^^gSk-Q + Ok- 



(22) 



Here Ok describes non-condensed magnons with k ^ ±Q. 
The leading contribution to the energy (per site) is ob- 
tained by neglecting the non-condensed particles alto- 
gether. Written in terms of condensate densities pi_2 = 
|V'±Qpj the energy reads 



— = -Sfi{pi + P2) + ^^i{p\ 



pI)+T2PiP2- (23) 



The coefficients ri^2 will be defined momentarily below, 
but first it is useful to discuss the important physics 
described by this simple equation in some generality. 
When p < 0, the energy is minimized by a vacuum state 
with no condensate, pi.2 = 0. Increasing p (by low- 
ering h) to positive values leads to a condensed state. 
A double-Q condensate emerges for Fi > F2, when 
Pi ^ P2 — Sp,/{Ti -I- F2). The physical meaning of this 
is found by expressing the spin operators in terms of bo- 
son condensates. Using ^l'±Q = ^pie*^^-^ and writing 



± 0_ , we find 



S+ = 2S 



S 



2/i 



riH-F2 



Ti+Fa 



e*''+cos[Q-r + 6'_], 
cos^[Q-r + 6l_]. (24) 



Hence, the double-Q condensate describes a coplanar 
magnetically ordered state. Observe that the 5*^ spin 



nar state described by Eq. (24) thus represents a 'super- 
solid' phase of the magnet. Xsingle-Q condensate arises 
when Fi < F2; in this case pi — Sfx/Ti and p2 — 
or vice versa. This state corresponds to a conventional 
umbrella (or cone) magnetic structure: 




S+ = S•,/^cxp[iQ•r + ^6'l], 
1 1 



O™ 



(25) 



The classical (order l/S*") expression for Fi 2 is given 
by 

Fi = J(0) - J(Q), F2 = J(0) - 2J(Q) + J(2Q). (26) 

The accidental degeneracy of the classical isotropic tri- 
angular antiferromagnet discussed in Sec. |III A] manifests 
itself via the relation Fi - F2 = J(Q) - J(2Q) = for 
Q = (47r/3,0). Thus the coplanar and cone states are 
degenerate classically. At first order in 1/5* one finds 



ri>F2 



(27) 



so that quantum fluctuations select coplanar state, con- 
sistent with our earlier spin- wave analysis. The corre- 
sponding calculation is sketched in Refs. I45I48I and will 
not be discussed here for brevity. Instead, here we take 
Eq. ( 26 1 as given and ask how the addition of weak DM 
interaction in Eq. ([2]) as well as spatial anisotropy influ- 
ences the delicate balance between coplanar and umbrella 
states at high magnetic fields. 



B. Spatially Isotropic Model with 
Dzyaloshinskii-Moriya Interactions: 

J = J',D/0,h = /iz 

In the geometry with h — hz the DM Hamiltonian in 
Eq. ([2]) reduces to a simple quadratic form of magnon 
operators. 



iD 
~2 



Hdm - ^ Z^ [^r S^+g. - h.c. 

rj = l,3 

= iDS J2 (4ar+6j - h.c") (28) 

= -2NDS[am{q ■ di) + sin(Q • (53)](pi - P2). 

As expected physically, DM coupling breaks the symme- 
try between ±Q points, and lowers the energy of the 
ip+Q condensate (for positive D > 0). The ground state 
follows from minimizing 



E 

N 



-Sp{pi + P2) - Sg{pi - P2) + 

-T^lipl + pI) +^2PlP2, 



(29) 
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where we introduced g = 2\/3D > 0. Following Ref. | 
we parameterize the densities as 



Pi 

P2 



poos 

• 2 

psm 



and minimize the energy with respect to </> and p. 
Three possible solutions exist: 



(ii) '^= 2'^ 
(in) cos^ (f) - 

E _ 

N ~ 



Sjf^ + g) E^ 

Ti ' N 
_S(ji-gl E 

' N 



2ri 



Ti 



2ri 



1 



.g(ri + r2) " 
Ai(ri-r2) 



2fiS 



ri + r2 



(30) 

(31) 

(32) 

(33) 



Solutions (i) and (ii) describe umbrella states formed 
from a single-Q condensate. For positive g the umbrella 
corresponding to (i) with pi ^ 0, p2 = has lower en- 
ergy. The 'mixed' solution (iii), which only exists when 
-1 < g{Ti + T2)/p{Ti-T2) < 1, yields the lowest energy 
provided that Eq. ( p7| is satisfied. This solution repre- 
sents the distorted V state of Fig. [2jd), discussed in Sec. 
|IV A3| above. The condensates here satisfy 






pS 

ri + r2 

pS 



gs 
Ti - r2 ' 

.9^ 



Fi -|- F2 Fi — F2 
and, as a result, the spin structure is non-coplanar, 

S+ = V2^e"'+[(|V'+Q| + |V-Ql)cos(Q-r + 0_ 

+i(|7/.+Q|-|V-Q|)sin(Q-r + (?_ 



(34) 



(35) 



The degree of non-coplanarity is controlled by the DM 
coupling g. The distorted V state has finite longitudi- 
nal chirality [recall Eq. (Is])] which is proportional to the 
density imbalance between the ±Q condensates. 



«,~25(|V+q|'-|V'-q|'). 



(36) 



Thus, when the magnetic field orients along the DM 
vector, the transition from the fully polarized state pro- 
ceeds in two steps as summarized in the schematic phase 
diagram of Fig. [9] Initially, at p = —g, one enters the 
umbrella state with ordering wavevector Q. This state 
persists at lower fields until p = g{Ti + F2)/(Fi — F2), 
where it is replaced by a distorted V state with conden- 
sates at both ±Q wavevectors. These findings agree fully 
with a very different analysis in Sec. |IV A] which was car- 
ried out for specific values of D/J and 1/S. Our results 
here show that the two phases obtained at high fields 
near saturation in fact emerge very generally provided 
DM coupling is weak and quantum fluctuations can be 
adequately captured perturbatively. 
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FIG. 9: (Color online) Phase diagram of the model in Eq. 
(29 1. The fully polarized state corresponds to p — 0, while all 
other states exhibit a finite boson density, p > 0. Umbrella 
states correspond to 4> = and 7r/2, while the mixed state 
with < (j> < tt/2 represents a distorted V phase. Observe 
that the origin of the phase diagram represents a tetracritical 
point. 



C. Spatially Isotropic Model with 
Dzyaloshinskii-Moriya Interactions: 

J = J',D/0,h = ftx 

To apply the formalism developed in this section to 
the case where the field orients perpendicular to the DM 
vector, it is simplest to rotate D and h such that h = /iz 
and D := Die. The DM Hamiltonian then reduces to the 
following combination of three magnon operators. 



- alar{ar+s, ^4+6,) 
Focusing on the condensates this expression reduces to 

/2S 



(37) 



^DM = D 



2i 

-i2QSj 



E 



\i}*+Q^P-Q{i}- 



'/'+Q)e" 



-i3Qr 



e-''J"'^] + {Q^-Q} 



(38) 



At this point, commensurability of the three-sublattice 
spin structure acquires crucial importance, since this im- 
plies e**'^'^'' — e^"^^^^ = 1 for all sites r of the triangular 
lattice. The end result, obtained by separating the mag- 
nitudes and phases of the condensates using Eqs. (30), is 
rather compact: 



i^DM = N^\f%SDp^l'^ sin 61+ sin(36'_). 



(39) 



Notice that Eq. (39 



is independent of (/>. This angle 
is then determined so ely by quantum effects which se- 
lect coplanar order, the detailed structure of which is set 
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by the angles 9±. In principle, quantum fluctuations do 
distinguish the latter angles and favor V ordering near 
saturation; this selection is, however, very weak and ap- 
pears only when one includes higher-order terms in Eq. 
(23) (involving six boson fields)^^. In contrast DM cou- 
pling distinguishes these angles more readily. One can 
immediately observe that, for D > 0, minimization of 
Eq. (39) requires that 9+ = ±7r/2. Let us choose for con- 
creteness 9+ = +7r/2. By consulting Eq. (24), one finds 



that the plane in which spins order is orthogonal to the 
(rotated) DM vector. That is. 



5? = 

Sy = ^/8Spcos{Q■r + 9^). 



(40) 



At the same time, we need to impose sin(30_) = —1 
given our choice for 9^. The different solutions for 6'_ 
simply describe equivalent magnetic structures connected 
by lattice translations, so we focus for concreteness on 
9_ = tt/2. This choice results in S^ = -y/8Spsm{Q ■ r). 
For the triangular lattice the product 



Q 



2tt 



-V mod 27r 



(41) 



takes on three inequivalent values: (y = 0), 27r/3 
{ly = 1), and — 27r/3 {v — —1). Correspondingly, the or- 
dered S^ spin components take on values which are zero 
{f = 0), positive (i/ = 1) and negative {i^ — —1). This 
represents the three-sublattice inverted-Y state found 
previously in Section IV B for a specific value of D/J 
and 1/S. The rather different approach adopted here re- 
veals that the onset of inverted Y order induced by DM 
coupling near saturation is in fact a very generic conclu- 
sion. 



D. Interplay Betw^een Spatial Anisotropy and 
Dzyaloshinskii-Moriya Interactions: J ^ J' ,D ^ Q 

We will now begin addressing for the first time the in- 
fluence of spatial anisotropy, which can be incorporated 
in a particularly simple manner near the saturation field. 
Let us first address the classical model without DM cou- 
pling. When J 7^ J' the Fourier transform of the ex- 
change interaction in Eq. ( 20 ) is now minimized at gen- 



erally incommensurate wavevectors ±Q given by 



Q = 2 cos 



-J' 
2J 



(42) 



Consequently the difference between Ti^2 defined in Eqs. 
(26) is non-zero already at the classical level; to order 
(J- J')VJ^ one finds Fi - Fs = -9( J - J'f/{2J). It 
follows from our analysis in Sec. |V A| that since F2 > Fi 
here, incommensurate umbrella states are stabilized clas- 
sically below the saturation field, even with arbitrarily 
weak exchange anisotropy. 

The inclusion of quantum fluctuations changes the sit- 
uation in an interesting way. Recalling that Fi > F2 due 



quantum effects in the isotropic limit, to leading order in 
1/S and second order in exchange anisotropy we have 



Fi 



aj 



9( J - J' 
2J" 



/\2 



(43) 



for some constant a. It is now apparent that the planar 
ordering favored by quantum fluctuations is in fact sta- 
ble over a finite range of anisotropy, though due to the 
shift in Q such planar configurations will now be incom- 
mensurate. Only when the anisotropy reaches a critical 
strength are these states supplanted by umbrella order. 

This analysis can be readily extended to incorporate 
DM coupling as well. Consider first the field orientation 
h = hz, where the magnetic field and D vector are par- 
allel. As long as spatial anisotropy is sufficiently weak 
that Fi exceeds F2, the phase diagram shown in Fig. 
B^ remains qualitatively intact (though the umbrella and 
distorted V states become incommensurate) . As \J — J'\ 
increases, the incommensurate distorted V state arises 
over a progressively smaller region of the phase diagram 
until, when Fi = F2, it is removed entirely. At larger 
anisotropy only incommensurate umbrella phases appear 
just below saturation. 

The case where the magnetic field orients perpendicu- 
lar to the DM vector is even more straightforward. In the 
preceding subsection the effectiveness of DM coupling re- 
lied critically on the wavevector Q being commensurate; 
see Eq. (38). With incommensurate Q the correspond- 
ing term sums to zero and thus drops out. Thus the 
DM interaction is effectively gone and we then recover 
the physics discussed above for the anisotropic quantum 
model without DM interactions. 



VI. SPATIALLY ANISOTROPIC MODEL: 

A. Classical limit 

In the remaining sections we will endeavor to address 
the phase diagram in the presence of spatial anisotropy 
at arbitrary fields. We start by treating the case with- 
out DM interactions for simplicity. Arbitrarily weak 
anisotropy J— J' destroys the accidental classical ground 
state degeneracy of the isotropic model in favor of an in- 
commensurate spiral ground state at all fields, similarly 
to the situation described above near saturation. The 



ordering wavevector is in fact given by Eq. ( 42 ) for all 



values of the magnetic field below saturation; the corre- 
sponding spin state is described by 

S(r) == sin 9 [cos (Q • r)x + sin (Q • r)y] -I- cos 9z. (44) 

While a unique classical ground state therefore emerges 
for all fields at zero temperature, the finite-temperature 
phase diagram is much more complicated and interesting. 
This can be anticipated on physical grounds. Indeed, as 
a general rule thermal fluctuations prefer coplanar over 
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FIG. 10: (Color online) Phase diagram of the classical 
anisotropic model with J' j J — 0.765. Interestingly, at this 
anisotropy strength the entropically favored planar phases are 
absent (within the resolution of our numerics), save for the 
coUinear UUD state which now requires intermediate temper- 
atures to appear. 



non-coplanar states (see the discussion in Sec. Ill); thus 
for sufficiently weak anisotropy one can expect that the 
entropic gain from planar order should be able to over- 
come its classical energy cost over a range of tempera- 
tures. This means that planar states should appear above 
some critical temperature Tpi. It is reasonable to expect 
that Tpi ~ {J — J'Y / J since the classical energy gain by 
umbrella states occurs at second order in J — J'. It is 
also clear that as J — J' increases, the planar states will 
be gradually pushed towards higher temperatures and 
disappear altogether above some critical value (which is 
magnetic-field dependent) of the spatial anisotropy. 

Of all the planar phases considered the UUD state, 
which in fact is coUinear, yields the highest entropy at 
finite temperature. The relative stability of the UUD 
state over the Y and V states is clear already at the spa- 
tially isotropic point: Figure [3] shows that at finite T 
the coUinear state expands at the 'expense' of the pla- 
nar ones. We thus expect the UUD state to persist the 
most upon deformation of the model's parameters. Our 
numerical findings fully support this conclusion and re- 
veal the rather non-trivial phase structure of the classical 



model. Figure 10 displays the rough phase diagram ob- 
tained using classical Monte Carlo with J' / J = 0.765. 
This value was chosen due to its closeness to the esti- 
mate for Cs2CuBr4, and the fact that for this anisotropy 
strength the wavevector Q 'fits' into the 48 x 48 system 
so that incommensurate orders are not frustrated in the 
geometries we simulated. In these simulations we have 
taken 20, 000 Monte Carlo steps for thermalization and 
20, 000 more for measurements. 

Remarkably, the UUD state indeed remains and 'floats' 
above the energetically preferred umbrella phase. The 
UUD order can be clearly identified by the behavior of 



dM/dh versus field. As shown in Fig. 11 at an intermedi 



FIG. 11: (Color online) dM/dh versus h aX T ^ 0.168J 
(green circles) and T = 0.12J (blue squares) for the classi- 
cal anisotropic model with J' / J = 0.765. 



ate temperature of T = 0.168 J we observe a pronounced 
dip in dM/dh, indicating a diminished slope of the mag- 
netization in the field interval between hd and hc2- Of 
course slower growth of the magnetization is expected 
for the magnetization plateau. At a lower temperature 
of T = 0.12J, dM/dh shows no such variation indicat- 
ing the absence of UUD order. This picture is further 
corroborated by the temperature dependence of the chi- 
rality as Fig. [12] illustrates for h = 2. 3 J. One clearly sees 
a discontinuous jump at T « 0.13 J in both the trans- 
verse and longitudinal chiralities, indicating a first-order 
finite-temperature transition from the entropy-stabilized 
UUD state to an umbrella phase. Note that, with the 
possible exception of the small regions near the UUD 
boundary where our numerics do not have enough accu- 
racy to reach any definite conclusions, the Y and V planar 
states are absent in the phase diagram. At the anisotropy 
we analyzed here they are replaced by the energetically 
favorable umbrella structure. 

The already non-trivial phase diagram we obtained 
here certainly deserves more extensive numerical in- 
vestigation. In particular it would be interesting to 
perform a systematic study increasing J — J' from 
zero to explore the collapse of the Y/UUD/V states 
and accompanying onset of umbrella order. We note 
that a similar finite-temperature competition between 
energetically-favorable and entropically-favorable states 
has been previously report ed in more complex spin sys- 
tems: frustrated pyrochloreSSEB and Shastry-Sutherland 
antiferromagneta^. It is worth noting that the roots of 
this behavior can be traced to the famous "Pomeranchuk 
effect" in '^He where the crystal phase of "^He, which is 
characterized by exponentially weak in distance exchange 
interaction between localized spins, has higher entropy 
than the normal Fermi- liquid phase. As a resu lt, upon 
heating the liquid phase freezes into a solicpSIMl 
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FIG. 12: (Color online) Chiralities k (blue squares), k_l (red 
triangles) and k^ (green circles) versus T at ft = 2.3J for the 
classical anisotropic model with J' j J — 0.765. The discon- 
tinuous jumps near T = 0.13 J signify a first-order transition 
between umbrella and UUD phases. 



B. Pseudo-quantum Ground States 

Our goal now will be to use Monte Carlo numerics 
to understand the zero-temperature phase diagram when 
both exchange anisotropy and quantum effects (modeled 
via the biquadratic approximation discussed in Sec. Ill E ) 



are present. Again to avoid becoming trapped in a local 
energy minimum, the phase diagram was obtained using 
simulated annealing as described in Sec. |IIIE} 

To explore the competition between anisotropy and 
quantum effects, simulations were performed with 
fixed S = Xjl but numerous anisotropics J' j J = 
1,0.95,0.9,0.85,0.8,0.765,0.7. We used M,dM/dh,K 
and K to construct the phase diagram at each J' /J; Fig. 
13 summarizes our findings. An interesting feature of the 
phase diagram is the complete absence of non-coplanar 
umbrella states, despite the fact that this type of or- 
der uniquely minimizes the classical energy at all fields. 
Even for rather substantial anisotropy J — J' ^ 0.3J, 
quantum effects (modeled here via the biquadratic inter- 
action) qualitatively alter the magnetization process. 

A second interesting feature is the persistence of 
commensurate Y/UUD/V states despite the exchange 
anisotropy. The UUD state is particularly stable, and at 
least within the approximations used is hardly affected by 
the finite J — J' values studied. In its vicinity commen- 
surate Y and V states appear over a finite field interval, 
with the latter being more robust than the former against 
anisotropy. All of this is in agreement with the previous 
large-S* analytical investigation of Ref. |23l Incommensu- 
rate planar phases, which reflect a nontrivial compromise 
between quantum fluctuations and anisotropy, appear at 
low flelds and near saturation; in the thermodynamic 
limit, these are expected to occupy a progressively larger 
portion of the phase diagram as anisotropy increases. 

It is worth briefly remarking on finite-size effects in 
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FIG. 13: (Color online) Ground state phase diagram for the 
spatially anisotropic triangular antiferromagnet with quan- 
tum fluctuations modeled via the biquadratic approximation 
discussed in Sec. |IIIE| All data points were obtained with 
Monte Carlo simulated annealing numerics. Dashed lines rep- 
resent phase boundaries which are absent due to finite-size ef- 
fects in our simulations but are expected on general grounds 
for an infinite system. 



our simulations. First, in Fig. 13 the high-field incom- 
mensurate phase appears over a broader field range at 
J' I J = 0.765 than J' j J = 0.7. This artifact arises be- 
cause the incommensurate spin structure that would ap- 
pear in an infinite system is less frustrated by the periodic 
boundary conditions at the former anisotropy strength. 
Second, we note that incommensurate order sets in at low 
and high fields only when J' < 0.85 J in our simulations. 
Closer to the isotropic limit, they are simply absent. This 
reflects finite-size effects arising from the relatively small 
systems modeled here. We expect that the phase bound- 
ary between the Y and incommensurate planar states em- 
anates from the h = 0, J' = J point. Similarly, the phase 
boundary between the V and the incommensurate planar 
states extends all the way to h — hse,t,J' — J', indeed. 



in Sec. VD we found that arbitrarily weak anisotropy is 
sufficient to produce incommensurate planar order near 
saturation. These phase lines represent commensurate- 
incommensurate transitions. We have indicated the ex- 



pected transitions with dashed lines in Fig. 13 



VII. SPATIALLY ANISOTROPIC MODEL WITH 
DZYALOSHINSKII-MORIYA INTERACTIONS: 

J/ J',D/0 

A. Classical Ground States 

Finally, we are in position to consider the full model 
featuring all of the ingredients we set out to study: spa- 
tial anisotropy, DM coupling, and quantum fluctuations. 
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FIG. 14: (Color online) Classical phase diagram for spins with 
spatially anisotropic exchange interactions, DM coupling, and 
a magnetic field applied perpendicular to the DM vector. 
Phase boundaries were computed analytically to leading order 
in DM strength D/J and anisotropy strength (J — J')^ /J^ as 
described in Appendix [B] 



As a first step we will establish the phase diagram in the 
classical case. For fields h = hz directed along the DM 
vector the classical ground states simply correspond to 
incommensurate umbrella order for all fields up to sat- 
uration. This outcome is extremely natural given our 
earlier findings that for this field orientation both DM 
coupling and anisotropy separately favor umbrella states. 



The phase diagram is more subtle for fields h = /ix 
oriented perpendicular to the DM vector. Wc found ear- 
lier that spatial anisotropy favors incommensurate um- 
brella states, while DM interactions prefer inverted Y or- 
der; thus, the resolution of their competition is far from 
obvious. Focusing on J' /J — 0.765 (again, this value 
minimizes finite-size effects) and D/J = 0.05, we find 
using simulated annealing that incommensurate copla- 
nar order arises at low fields, h < 0.34/isat, reflecting 
a nontrivial compromise between these competing inter- 
actions. DM coupling dominates at intermediate fields 
0.34/isat '^ h < 0.8/isati where inverted-Y states appear. 
The appearance of such a broad commensurate state in 
the anisotropic system even with only quadratic spin cou- 
plings is rather remarkable. Finally, at larger fields up to 
saturation spatial anisotropy dominates, leading to non- 
coplanar umbrella order. One can in fact analytically 
estimate the phase boundaries between these three spin 
states found in our numerics. The calculation is described 
in Appendix [B] and the resulting classical phase diagram 
appears in Fig. [M] 



B. Pseudo-quantum Ground States 

Let us now explore how quantum effects — again mod- 
eled within the biquadratic approxiation — modify the 
classical phase diagrams discussed above. Consider first 
the h. = hz field orientation. The situation at high fields 
was already analyzed in Sec. |VD[ where we found the 
presence of umbrella order just below saturation, followed 
by an incommensurate planar state provided anisotropy 
was not too strong. At intermediate fields we found in 
the quantum problem with either DM coupling or spatial 
anisotropy that commensurate planar phases emerged. It 
is thus reasonable to anticipate the same outcome when 
both elements are present, at least for sufficiently weak 
DM coupling and anisotropy. At low fields DM cou- 
pling led to umbrella order in the quantum problem an- 



alyzed in Sec. IV A3 while the interplay between spatial 
anisotropy and quantum effects led to an incommensu- 
rate planar state in Sec. |VIB[ 

Putting together these findings suggests that the phase 
diagram for the spatially anisotropic pseudo-quantum 
model depicted in Fig. [13] evolves in the following man- 
ner as one increases the DM coupling strength from zero. 
First, umbrella order immediately begins to 'eat away' at 
the planar phases just below saturation, occupying a pro- 
gressively larger fraction of the high-field phase diagram 
as the DM coupling increases. The low-field incommen- 
surate planar state stabilized by the interplay between 
spatial anisotropy and quantum fluctuations is more ro- 
bust against DM coupling. Only beyond a critical value 
of the DM coupling does umbrella order begin to take 
over in the low-field portion of the phase diagram. Of 
course more complicated scenarios are all possible, par- 
ticularly if DM coupling and/or anisotropy are not espe- 
cially weak; a detailed study of the problem for this field 
orientation would be interesting to carry out in future 
work. 

Our main focus, however, is on the low-symmetry field 
orientation h — /ix where the magnetic field and DM 
vectors are orthogonal. This is the physical situation rel- 
evant for the interesting experiments of Ref. |T9l which 
motivated this study. We explored the zero-temperature 
phase diagram here using extensive simulating annealing 
numerics, modeling quantum fluctuations as before using 
the biquadratic approximation. We note that particular 
care must be taken when performing these simulations 
to avoid spurious finite-size effects. In particular, at low 
fields in 48 x 48 systems we found an unusual incom- 
mensurate planar state exhibiting structure-factor peaks 
at two incommensurate wavevectors with non-zero mo- 
mentum Qy along the y-direction. This phase, however, 
proved to arise due to finite-size effects — upon increasing 
the system size to 192 x 192, order characterized by a 
single wavevector and vanishing Qy emerged. 

Figure 16 summarizes our results. As noted in Sec. 



|IIIB[ the Hamiltonian no longer possesses any contin- 
uous symmetries; hence some of the phase boundaries 
discussed earlier disappear from the phase diagram and 
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FIG. 15: (Color online) Magnetization versus field at tem- 
perature T — O.OOIJ for the spatially anisotropic antiferro- 
magnet with J' — 0.765J, DM interactions satisfying D _L h, 
and quantum effects modeled via a biquadratic interaction. 
The curves shown correspond to different values of the DM 
strength: D = (blue circles), D — O.OIJ (green trian- 
gles), D = 0.02 J (red stars), D = 0.04 J (cyan squares), and 
D = 0.05 J (pink hexagons). The quick rounding of the lower 
edge of the plateau reflects the symmetry equivalence of the 
Y and UUD states when DM interactions are present in this 
fleld orientation. 



become crossovers. One of these is the transition be- 
tween the Y and the UUD states. The difference be- 
tween these phases in the problem without DM coupling 
originates from the finite superfluid component of the Y 
state, present due to spontaneous breaking of U(l) spin 
rotations about the field axis, along with rotation sym- 
metry which the Y state breaks but the UUD state does 
not. In the D _L h problem these two states are symme- 
try equivalent as they only break the same discrete lattice 
symmetries, and are thus not distinct phases. This shows 
up in our numerical simulations as a quick rounding of 
the lower end of the magnetization plateau upon increas- 
ing the DM coupling at fixed J' /J, as illustrated in Fig. 
[T5| Simultaneously, for D as small a.s D = O.OIJ the 
coplanarity K becomes finite, although small, inside the 
former UUD interval between the Y and V states. 

We also observe persistence of the distorted V state at 
intermediate fields above the (former) UUD state. [Note 
that as discussed in Sec. |IIIB| the symmetry distinc- 
tion between the distorted V and UUD orders persists 
in this field orientation, so that a bona fide phase tran- 
sition separates these states.] In fact it appears that the 
distorted V state is the most stable of all the commen- 
surate states considered previously. Because the spins 
in this phase can smoothly adjust to gain DM energy, 
this state survives even at the strongest DM coupling 
strength D/J = 0.05 considered by us. By contrast, 
the previously robust UUD state, having lost its symme- 
try distinction from the less stable Y state as discussed 
above, is seen in Fig. [15] to essentially disappear for such 
a strong DM coupling. 
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FIG. 16: Sketch of the phase diagram for the quantum spa- 
tially anisotropic model with varying DM coupling strength, 
J' /J — 0.765, and a magnetic fleld h oriented perpendicular 
to the DM vector. Data points indicate phase boundaries de- 
termined using simulating annealing numerics, with quantum 
fluctuations modeled via a biquadratic interaction. Dashed 
lines interpolate between these data points and are drawn for 
convenience. The location of the inverted Y phase is approx- 
imate as discussed in the main text. 



Our simulations also reveal that for sufficiently strong 
DM coupling D > 0.04J a narrow region of the com- 
mensurate inverted Y state appears above the distorted 
V phase. This is quite consistent with the phase dia- 



gram of Fig. 14 for the classical model described in the 
previous subsection: being a prominent phase there, the 
inverted Y state is also natural in the quantum prob- 
lem once quantum effects are sufficiently "weakened" by 
spatial anisotropy and DM interactions. The transition 
between the inverted Y and distorted V states appears 
continuous in our numerics and hard to pin down pre- 
cisely, in part because the overall extent of this phase 
is rather narrow. For this reason its phase boundary in 
Figure [16] is less accurate than for the other phases. 

The remaining high- and low-field regions are found 
to be occupied by the incommensurate coplanar states 
which owe their stability to quantum fluctuations. This is 
particularly clear for the high-field region near saturation 
where an incommensurate analog of the V state wins over 
the classical umbrella state only due to 1/5 interactions 
between spin waves, as discussed in Sec. |VD[ 



VIII. CONCLUSIONS 

We have explored the phase diagram of a spatially 
anisotropic triangular lattice quantum antiferromagnet 
subject to asymmetric DM interactions and an external 
magnetic field. By treating spatial anisotropy (J— J')/ J, 
quantum fluctuations due to the flnite spin value S, and 
DM coupling D as perturbations of the well-understood 
isotropic classical antiferromagnet, we have found a rich 
variety of behaviors sensitive to the relative strengths of 
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the perturbations considered. The root of this richness 
hes in the large accidental degeneracy of the unperturbed 
model with J' — J, S — oo and D ~ 0, along with the 
fact that each of these perturbations favors different or- 
dered states in the manifold of accidentally degenerate 
configurations. 

Our main findings are as follows: 

1) In agreement with numerous previous studies, for 
the isotropic model without DM coupling we observe 
that quantum fluctuations select coplanar Y and V states 
and a coUinear UUD phase out of infinitely many de- 
generate states available in the classical limit 5 — > oo. 
Quantum effects most effectively split this degeneracy 
in the vicinity of one-third magnetization where nearly 
coUinear low-energy states are accessible. This implies 
a greater robustness of the quantum-selected states to 
additional perturbations in these regions of the phase di- 
agram, a trend that is indeed seen throughout our study. 
We also established that quantum effects can be semi- 
quantitatively modeled within a purely classical Hamilto- 
nian by incorporating a biquadratic spin interaction with 
a field-dependent coupling chosen to mimic l/S correc- 
tions to the classical ground state energies. Such a purely 
classical model has the virtue of allowing conventional 
Monte Carlo simulations to be employed to ascertain the 
(approximate) phase diagram when all the competing in- 
teractions of interest are present. 

2) Adding DM interactions introduces two new com- 
mensurate states into consideration — distorted V and in- 
verted Y orders. With fields applied perpendicular to the 
triangular lattice plane, DM coupling stabilizes umbrella 
order classically at all fields up to saturation. Quan- 
tum effects, however, stabilize a distorted V state at in- 
termediate fields, bordered on both sides by the classi- 
cally driven umbrella phase. For all in-plane fields up 
to saturation, DM interactions select inverted Y states 
classically. Quantum fluctuations arc still more effective 
at modifying the phase diagram in this field orientation, 
producing Y order at low fields and a distorted V phase 
at intermediate fields. Only at high fields does the clas- 
sically favored inverted Y state appear in the quantum 
problem. It is worth pointing out here that the new dis- 
torted V and inverted Y states also appear naturally in 
the theoretically simple limit of high magnetic fields near 
saturation, as discussed in Sec. |V] 

3) Spatial anisotropy, on the other hand, prefers non- 
coplanar umbrella (cone) configurations classically. Even 
this simple classical model, which has a unique ground 
state, harbors some surprises. We find that for suf- 
ficiently small spatial anisotropy an entropic order-by- 
disorder selection prevails over energetic considerations 
over a range of temperatures and magnetic fields. This 
results in an abrupt first-order phase transition from an 
umbrella state into a fluctuation-stabilized UUD state at 



finite temperature; see Sec. VIA This physically reason- 
able but so far unexplored feature of the classical Heisen- 
berg model with anisotropic exchange interactions de- 
serves more extensive numerical investigation on its own. 



4) The competition between spatial anisotropy and 
quantum effects (simulated within an effective classical 
model with biquadratic spin couplings) results in a rich 
phase diagram shown in Fig. |13[ These competing ingre- 
dients compromise to form incommensurate versions of 
the Y and V states at low and high fields, respectively. 
Remarkably, the intermediate-field region of the phase di- 
agram in Fig.flSlis only weakly affected by the anisotropy 
strengths we analyzed, with the commensurate Y, UUD, 
and V states favored by quantum effects all appearing 
prominently. In particular, the UUD state which under- 
lies the one-third-magnctization plateau, shows no ob- 
vious reduction in its width. All these features of the 
phase diagram agree well with previous analyticaP^l and 



numericaP^ studies. 

5) When spatial anisotropy, DM interactions, and 
quantum effects (again incorporated in an effective clas- 
sical Hamiltonian) are all present, we obtain the phase 
diagram of Fig. |16[ which contains most of the phases re- 
viewed above. It is interesting to compare this with the 
experimental phase diagram of Ref. 19" which features up 
to nine different phases. Clearly our phase diagram is 
less diverse but still contains, for D « 0.04J, five differ- 
ent phases, most of which are coplanar. In particular it is 
tempting to speculate that one of the puzzling high-field 
phases (B, III, IV or 2/3, in the notation of Ref. ITO|) can 
be identified with the inverted Y state appearing near 
the phase boundary between the distorted V and the in- 
commensurate planar phases. Further experiments that 
directly probe the spin structure in the various phases 
seen are, however, required to make a more definitive 
comparison with our results. 

We emphasize here that nowhere in our study did we 
find hints of a possible two-thirds-magne tization plateau 
suggested by several experimental paperai^ElMl, in this 
regard the possibility of this novel magnetization plateau, 
at present identified theoretically only once^, requires 
more careful investigations. It must also be kept in mind 
that our treatment of quantum fluctuations is only per- 
turbative, and certainly in the S — 1/2 model at least 
the phase boundaries will be quantitatively different from 
what we have found. Whether or not strong quantum 
fluctuations can induce entirely new states, unseen in 
our semiclassical analysis, remains an interesting open 
question. 

More broadly, our analysis relies on the observation 
that near the spatially isotropic point, where the rele- 
vant low-energy states form a three-sublattice pattern, 
the only relevant DM interaction is that given in Eq. 
([2]). This implies that the phase diagram should be in- 
sensitive to the orientation of the magnetic field in the 
b — c plane of the material. Recent experimental find- 
ings support this statemenlPS, and, in our view, should 
be interpreted as lending strong support to our perturba- 
tive approach to Cs2CuBr4. It is worth pointing out that 
this is certainly not the case for the isostructural material 
CS2CUCI4; there it has been established experimentally 
that phase diagrams for the cases h || b and h || c are 
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indeed difFerentfl^. That difference has been attributed 
to additional symmetry-allowed DM interactionalSI and, 
fundamentally, is tied to the fact that CS2CUCI4 is very 
much a quasi-one-dimensional material, in contrast to 
Cs2CuBr4. 

Another, and so far little explored, way to gain insight 
into the rich physics of frustrated antiferromagnets is pro- 
vided by a magnetic system's response to various impuri- 
ties. Cs2CuBr4 has been probed in this way^H, and initial 
theoretical investigations have appeared recently!^ 

We hope that our study will stimulate further theo- 
retical and experimental studies of Cs2CuBr4, as well as 
related materials and models, in the near future. 
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Appendix A: Symmetry-allow^ed 
Dzyaloshinskii-Moriya interactions in Cs2CuBr4 

The general form of the DM interaction, consistent 
with crystal symmetry of Cs2CuBr4, was derived in Ref. 
[T51 for the isostructural material CS2CUCI4. For the 
three-dimensional crystal with layers indexed by z, it 
reads 

HbM = ^(D' • Sr X Sr+5, + D+ • Sr X S^+S, 
r 

- D- • Sr X Sr-63) (Al) 

where 

D' = 2D,{-irz + Dy{-l)yy, (A2) 

D± = D',{~irz±DU~ir+'& + D'yi-l)y+^yiA3) 

We are, however, interested in a 2D triangular antiferro- 
magnet so will henceforth consider only one layer with 
z — for concreteness. We now show that for a weakly 
deformed triangular antiferromagnet, which is an appro- 
priate characterization of Cs2CuBr4, only two DM cou- 
plings Dz and D'^ need to be retained to leading order 
in perturbation theory. Furthermore, D^ and D'^ yield 



first-order contributions which are in fact identical in na- 
ture and can thus be combined into a single DM coupling 
oriented along the crystal z-axis. 

These results follow from the observation that all clas- 
sical ground states of the isotropic model form a three- 
sublattice structure which repeats in the y-direction; see 
Fig. fl] Let us first consider the term Dy. Because of 
the repeating structure in the y-direction, we see that 
the cross products of spins from the A, B and C sublat- 
tices cancel out due to the oscillating (—1)'' factor in the 
Dy term. Thus Dy can only appear at second-order or 
higher in perturbation theory, and for this reason can be 
neglected. Similar reasoning dictates that D'^ and Dy, 
which also include oscillating (—1)^ factors, do not con- 
tribute to first order and can thus be omitted as well. 

Consider now the remaining terms, D^ and D'^. From 
Fig. [1] we see that Sr+^j and Sr-Sj correspond to the 
same sublattice. Likewise, spins Sr+^a ^^^ ^r-<53 can 
also be identified. Consequently, the most general DM 
Hamiltonian above, when evaluated in an arbitrary three- 
sublattice classical ground state, reduces to the following 
expression, 

r 

- -D^SrXSr-Sa) 

= ^ Z • [{D,Sr X Sr+6, - D'^S^ X S.^Ss) 
r 

- {D,S, X Sr-s, - D',S, X S,+sJ] 

^ {D, ~D',)J2^- (Sr X Sr-53 - Sr X S,+s,) 
r 

= (i^; - Z?.) ^ Z • [Sr X (Sr+5, + Sr+,3)].(A4) 
r 

Hence, defining D — {D^ — D'^)i = Dz, we arrive at the 
DM Hamiltonian in Eq. ^. 



Appendix B: Analytical determination of phase 
boundaries in Figure |14| 

In this Appendix we will sketch the derivation of the 
phase boundaries depicted in Fig. [14] We adopt a 
variational approach and calculate the energies for the 
three phases found in simulations — incommensurate pla- 
nar, commensurate inverted-Y, and incommensurate um- 
brella states — to determine which minimizes the energy 
as a function of field. Of course there is no guarantee that 
only these three phases are relevant, but our numerical 
findings suggest that this is the case for weak anisotropy 
and DM couphng. 

For concreteness, let us take the field in the y-direction. 
For the incommensurate planar state, we work at low 
fields and parametrize the spins as 

Sr = cos[Q • r + 0(r)]x + sin[Q • r + 0(r)]y, (Bl) 

where Q is the incommensurate wavevector from Eq. 
(42 1. With 0(r) = 0, this expression yields the exact 
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classical ground state at zero field and in the absence of 
DM coupling. [We will neglect the dependence of Q on 
the DM coupling here, which results in energy correc- 
tions of order (D/J)^.] In finite fields the function 0(r) 
becomes non-zero and gives rise to a net magnetization 
along the y-direction; it suffices to take 



by Eq. ( 44 1 ; the corresponding energy per site is 

J'f 



N 



2 J' + J'- 
2J 



hM - 



(2 J 



2J 



M" 



(B9) 



(r) = 2M cos Q • r. 



(B2) 



Up to order h? , M represents the magnetization which is 
given by 



To leading order in D/J and (J— J')'^/J^, the incom- 
mensurate planar and inverted-Y state energies cross at 
a magnetization Mi satisfying 



M 



hJ^ 



2J 



(I-IOM^) = %/3D[(l - Mi)=^/Vl + 3Afi 



(2 J + J')[4J3 + J'3 _ 2J2(J' - D^4 - [J'/Jf 



(B3) 



(This expression can be obtained by evaluating the en- 
ergy to second order in /i, and then minimizing the ex- 
pression to find M .) The incommensurate planar energy 
per site, to order K^, is then 



E. 



planar 

N 



D 



hM + — v/4-(J7J)2[J - (2 J + f)M^] 
(2J + J')2(2J2 - 2JJ' + J'2)Af2 



2 J2 + J'^ 
2J 



2J3 



- l + 3Mf]. (BIO) 

For M < Ml the planar state has lower variational en- 
ergy, while for M > Mi the inverted-Y state wins. Tak- 
ing D/J = 0.05 and J' /J = 0.7, we find Mi « 0.3. Al- 
though our calculation of the planar energy was perturba- 
tive in the field and is thus most reliable near M = 0, this 
approximate result for Mi agrees remarkably well with 
the lower phase boundary determined numerically for 
these parameter values. The umbrella state and inverted- 
Y energies balance at a magnetization M2 which satisfies 

• (j-jy 



(B4) 



2J 



(1 + M2) = D^3{l + 3M2){1-M2). (Bll) 



Both the inverted-Y and umbrella states can be easily 
found at arbitrary fields. For the inverted-Y state, we 
have 



Sa = y 
Sb = 

Sc = 



(B5) 



1 - - (3Af - l)2x + - (3M - l)y (B6) 



l-i(3M-l)2x+i(3Af-l)y, (B7) 



In this case the inverted-Y state yields a lower energy 
for M < M2. With .r/J = 0.7 and D/J = 0.05 we ob- 
tain M2 ~ 0.7, also in excellent agreement with numer- 
ics. At a critical anisotropy strength (J — J'Y / J ~ 3D, 
the magnetization values Mi and M2 coincide, and for 
larger anisotropy the inverted Y state no longer appears 
in the phase diagram. We then have only two phases 
which compete — the incommensurate planar and um- 
brella states. To leading order in D/J and (J— J'Y / J"^, 
their energies cross at a magnetization 



which yields an energy per site 

-LVinvcrtcd— Y 



N 



Uj + 2J') + hM - hj + 2J')M^ 



M.^ 



9{J-J'y 
2V^JD 



-1/2 



(B12) 



+ 1)^3(1 + 3Af)(l-Af)3/2 



(B8) 



The spin configuration for the incommensurate umbrella 



For M > M3 the incommensurate umbrella phase 
emerges, whereas for M < M3 incommensurate planar 
order appears. Plotting the phase boundaries Mi^2,3 de- 
rived above leads to the classical phase diagram of Fig. 
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